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Specification 

background 

1. Technical Field 

The present invention relates to high fidelity audio systems and, more particularly, high 

fidelity digital audio systems that require sampling rate modification. 

2. Related Art 

There exist different sampling rate conversion algorithms to alter the sampling rate of a 
signal from one sampling rate, fj, to another sampling rate, f2. Typically, it is essential to change 
between two different sampling rates without degrading any of the useful information in the 
original digital signal. Several techniques are used traditionally. The conceptually simplest 
method involves converting the digital signal sampled at the rate fi into an analog signal using an 
accurate digital-to-analog converter (DAC). The resulting analog signal is lowpass filtered to 
remove digital images and then sampled at the rate f2 using an accurate analog-to-digital 
converter (ADC). While conceptually simple, the harmonics and noise distortion occurring in the 
conversion between the analog and digital domains degrades the quality of the resulting signal to 
a degree that renders the method unusable for high fidelity applications such as digital stereo 
audio. 

A conventional multi-rate technique has been popularly employed for digital sampling 
rate conversion for a fixed ratio M/N, where M and N are positive integers. This method is only 
effective when M/N is quite large or small; for M/N « 1, a task referred to as "near-unity 
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fractional sampling rate alteration", the method demands excessive computational complexity. 
For example, a task of primary focus for the present invention is the conversion of digital audio 
between the sampling rates of fi=44.1 kHz and f2=48 kHz, a task commonly referred to as 
"compact disk (CD) to digital audio tape (DAT) conversion", as shown in Figure 1 . Expressed as 

5 a ratio of integers, M/N = 160/147, or, M/N « 1.088. The input signal is first up-sampled by a 
factor 160, then passed through a sharp digital lowpass filter, and finally decimated by a factor 
147. Implementing a 160-fold interpolation for digital audio requires a computationally very 
complex digital algorithm and this method is thus undesirable for CD to DAT conversion. 

What is needed, therefore, is a method and an apparatus that addresses the problem of 

10 efficiently altering the sampling rate of a high fidelity digital signal, such as digital audio, 
between two sampling rates that are only slightly different. 



2 



PATENT APPLICATION 
Docket No. BP 3373 

Summary of the Invention 
The present invention provides an efficient method for near-unity sampling rate alteration 
in high performance applications, such as CD to DAT conversion. Specifically, the input digital 
signal is first interpolated (up-sampled) by a factor 8 and lowpass filtered to form an 
5 intermediate signal. The up-sampled signal is then decimated by a factor of 7.35 in the 
described embodiment of the invention. Circuitry that implements a clamped cubic spline 
interpolator (CCSI) algorithm is then employed to accurately interpolate the intermediate signal 
to points in-between adjacent samples of the intermediate signal as required by the 48 kHz 
output sampling rate. The CCSI is highly accurate due to highly accurate derivative estimates 
10 arrived at by repeated Richardson extrapolation. In the example CD to DAT converter covered 
in detail, fourth order Richardson extrapolation is employed. It is shown by this example that 
the proposed method yields the desired performance, is computationally efficient and requires 
little storage. 

More specifically, a sampling rate converter for altering an initial sampling rate to a 
15 desired sampling rate resulting in near-unity fractional sampling rate alteration includes an 
interpolator for receiving and up-sampling digital data representing a continuous waveform 
signal. The interpolator produces upsampled and filtered data at an upsampled rate relative to an 
initial digital data sampling rate. An accumulator that receives a clock at the upsampled rate 
determines approximate values of time, x, that relate to sampling points at the desired sampling 
20 rate and produces a trigger at the sampling points at the desired sampling rate. Thus, the 
circuitry produces a function value (signal magnitude) whenever the trigger is generated. 

A signal buffer receives the upsampled and filtered data at the upsampled rate and 
produces buffered samples that are N samples wide wherein N is a value determined by a 



3 



PATENT APPLICATION 
Docket No. BP 3373 

Richardson formulation for generating derivative estimates. A derivative estimation block 
receives the buffered samples and produces at least two derivative estimates based upon the 
buffered samples. A coefficient calculation block then produces two polynomial coefficient 
values based upon the at least two derivative estimates produced by the derivative estimation 

5 block and also based upon at least two select buffered samples produced by the signal buffer. 
Finally, a polynomial evaluation block for producing a function value representing an estimated 
sampling point value produces the function value (estimated signal magnitude) based upon the 
polynomial coefficient values and the approximate value of (x). The derivative estimation block, 
the coefficient calculation block and the polynomial evaluation block produce outputs based 

10 upon the trigger produced by the accumulator. The trigger represents a sampling point at the 
desired sampling rate. 

More generally, the above described circuitry, in the described embodiment, includes an 
interpolator that up samples the digital data eight times and a clamped cubic spline interpolator 
that decimates the upsampled data by a factor of 7.35. Accordingly, a hi-fidelity signal having a 
15 44. 1 kHz sampling rate may be converted to the same signal having a 48 kHz sampling rate in a 
real time basis for hi-fidelity applications. 
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Brief Description of the Drawings 
A better understanding of the present invention can be obtained when the following 
detailed description of the preferred embodiment is considered with the following drawings, in 
.which: 

5 Figure 1 is a functional block diagram of a sampling rate converter according to one 

embodiment of the present invention; 

Figures 2A and 2B are signal illustrations that illustrate sampling and upsampling of an 
example signal and one problem addressed by the present invention; 

Figure 3 shows a top level block diagram of the present invention for application in CD 
10 to DAT conversion; 

Figure 4 shows the two- fold interpolation filters implemented in poly-phase form; 
Figure 5 shows an example implementation of LPF 3 in poly-phase form; 
Figure 6 shows the top level block diagram of the CCSI block; the intermediate signal, f, 
with sampling rate 352.8 kHz enters the block and it generates output samples at a rate of 48 
15 kHz; 

Figure 7 is a functional block diagram that illustrates a clamped cubic spline interpolator 
formed according one embodiment of the present invention; 

Figure 8 is a block diagram of the "Signal Buffer" appropriate for calculating derivative 
estimates using fourth order Richardson extrapolation according to one embodiment of the 
20 present invention; 

Figure 9 is a functional block diagram of the derivative estimation block according one 
embodiment of the invention; 
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Figure 10 is a functional block diagram of the CCS coefficient calculation block that 
shows the details of the calculation of the CCS coefficients according one embodiment of the 
invention; 

Figure 1 1 shows the details of the CCS evaluation block of Figure 7; 
5 Figures 12A and 12B are flow charts of a method of the present invention; 

Figure 13 shows the magnitude response in dB of the sampling rate converter in the 
frequency range 0-30 kHz; 

Figure 14 shows a close-up of Figure 13; 
10 Figure 15 shows the typical ripple of the passband magnitude response; 

Figure 16 shows a typical simulation result for one embodiment of the present invention; 
Figure 17 shows a block diagram of the scheme employed to evaluate the SNDR and 

SFDR; 

Figure 18 shows simulated SNDR and SFDR vs. frequency of the sinusoidal input signal; 

15 and 

Figure 19 shows the simulated magnitude response of the sampling rate converter found 
as the ratio of output signal amplitude and input signal amplitude. 
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Detailed Description of the Drawings 

Figure 1 is a functional block diagram of a sampling rate converter according to one 
embodiment of the present invention. In the described embodiment, the sampling rate converter 
of Figure 1 is used to convert a high fidelity audio signal with a sampling rate of 44.1 kHz to the 
5 same high fidelity audio signal with a sampling rate of 48.0 kHz. It is understood, of course, that 
the inventive approach may be used for other near-unity sampling rate conversion applications. 
Generally, converting from a first to a second sampling rate defined by a ratio of two or more is a 
relatively easy task. For example, one may readily convert by integer amounts merely by 
upsampling and lowpass filtering at a higher rate in a digital processor. As stated before, 

10 however, converting between nearly identical sampling rates by integer interpolation followed 
by integer decimation is a more complex conversion.. Such complex computational approaches 
are satisfactory for some circumstances, but are not satisfactory for real-time high fidelity 
applications. In the described embodiment of Figure 1, sampling at a 44.1 kHz sampling rate 
yields a signal that ranges from 0 to 22.05 kHz having a mirror image that ranges from 22.05 to 

15 44. 1 kHz. Because hearing range for high fidelity applications is from 0 to 20 kHz, the sampling 
rate of 44.1 kHz is adequate for high fidelity applications. One problem, however, is that 
different protocol devices sample at different rates though approximately at the same rate. For 
example, compact disc players typically sample at 44.1 kHz, while other digital applications 
sample at, for example, 48.0 kHz. Accordingly, to convert from a CD to a DAT requires near- 

20 unity sampling rate conversion. 

Table 1 below gives additional typical performance requirement details of the sampling 
rate alteration process. Many of these requirements are discussed further in the subsequent 
sections of this disclosure, but the following is a brief summary of the requirements. Specifically, 
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the signal passband is the frequency range 0-20 kHz. The ripple (variation) of the sampling rate 
converter magnitude response over the signal passband is limited to 0.01 dB. The phase response 
over the signal passband is required to be linear. The transition band, i.e., the region between the 
upper passband frequency and the lower stopband frequency is the frequency range 20 kHz-24 
5 kHz. The required stopband attenuation is specified to be at least 110 dB. For signals in the 
specified passband, a minimum of 96 dB signal-to-noise-and-distortion ratio (SNDR) should be 
provided, with a minimum of 100 dB spurious-free dynamic range (SFDR). 



Input Sampling Rate 


44.1 kHz 


Output Sampling Rate 


48.0 kHz 


Passband 


0-20 kHz 


Passband Magnitude Ripple 


< 0.01 dB 


Phase Response 


Linear 


Transition Band 


20 kHz-24 kHz 


Stopband Attenuation 


> HOdB 


SNDR 


>96 dB 


SFDR 


> 100 dB 
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Table 1 - Typical specifications for the "CD to DAT" sampling rate converter 

Figures 2A and 2B are signal illustrations that show sampling and upsampling of an 
example signal and one problem addressed by the present invention. Referring to Figure 2A, a 
15 random signal 02 is shown for which samples are taken at a specified sampling rate. As may be 
seen, random signal 02 is shown to be sampled 5 times at a specified sampling rate. A signal 04, 
which is a copy of signal 02, is shown below signal 02 with sampling points for a slightly higher 
sampling rate. As may be seen, the samples at the desired sampling rate of the signal 04 precede 
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the samples at the original sampling rate by increasing amounts. If the desired sampling rate 
were merely an integer multiple of the original sampling rate, the problem would be straight- 
forward. Merely increasing the upsample rate followed by appropriate lowpass filtering would 
result in samples being taken at the desired sampling rate. However, as may be seen, for near- 
5 unity sampling rate conversions, an approach for converting from an original sampling rate to a 
desired sampling rate is not as simple.. 

Referring now to Figure 2B, it may be seen that merely upsampling the original sampling 
rate may not be adequate. As may be seen, two sampling points are shown on a signal 06. On a 
signal 08, which is a copy of signal 06, the original sampling rate has been increased by a factor 
10 of 8. The sampling point on signal 08, however, still resides between the last two sampling 
points at the upsampled sampling rate. Accordingly, determining a signal magnitude of signal 08 
at the desired sampling point resulting from the desired sampling rate is not readily obtained 
merely by increasing a sampling rate, even by a significant increase in rate such as suggested in 
Figure 2B. 

15 As will be described in greater detail below, however, the invention takes advantage of 

the fact that increasing a sampling rate by a factor, for example 8 times, creates sampling points 
sufficiently close to a sampling point of a desired sampling rate, as shown in Figure 2B for the 
signal 08, that using the inventive techniques disclosed herein, the actual signal magnitude for 
the signal shown as a copy of signal 04 may readily be obtained with a high degree of accuracy. 

20 Figure 3 shows a top level block diagram of the present invention for application in CD 

to DAT conversion. It consists of three 2-fold interpolation filters to increase the sampling rate 
to 44.1x2 3 kHz = 352.8 kHz. Each interpolation stage conceptually consists of a block that 
increases the sampling rate 2-fold (the block labeled "T2") and then lowpass filters the resulting 
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signal (the block labeled LPF) to eliminate images resulting from the increase in sampling rate. 
For the example sampling rate converter, LPFi is of order 1 19, LPF 2 is of order 17, and LPF 3 is 
of order 1 1 . The final output samples must be calculated by interpolating in-between the samples 
of the intermediate signal with sampling rate 352.8 kHz, as dictated by the desired 48 kHz output 
5 sampling rate. As shown in Figure 3, this task is performed by a block referred to as "CCSI", the 
clamped cubic spline interpolator. 

The clamped cubic spline interpolator that decimates the upsampled sampling rate by a 
factor of 7.35 to produce a 48 kHz sampling data stream includes circuitry for implementing a 
plurality of algorithms in hardware to achieve the desired result. As was described above in 

10 relation to Figure 2B, the clamped cubic spline interpolator is required to highly accurately 
determine a signal magnitude (or function) value at a sampling point defined by the decimated 
sampling rate. One approach for achieving the overall result that is provided by the architecture 
of Figure 3 is to upsample by a factor of 160 and to decimate a factor of 147. This approach is 
computationally complex and time-consuming and is not satisfactory for real-time high fidelity 

15 applications. Upsampling by a factor of 8, interpolating using the CCSI algorithm, and 
decimating by a factor 7.35, as shown in the architecture of Figure 3, satisfactorily achieves 
signal-to-noise and distortion ratio (SNDR) requirements. The problem, then, is decimating by 
a factor of 7.35, which is the function of the clamped cubic spline interpolator. As an additional 
point, it would be straight-forward to merely upsample the original signal by a factor of 8. One 

20 drawback with that approach, however, is that a filter would have to be designed to have a much 
sharper response to satisfactorily remove all images resulting from the upsampling. 
Accordingly, upsampling occurs in three stages, wherein the sampling rate is doubled each of 
three times, and is then filtered by a digital finite impulse response (FIR) filter. 
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The sampling rate converter of Figure 3 performs a Taylor expansion and a Richardson 
extrapolation to get highly accurate derivative values for the two upsampled sampling points that 
appear about the sampling point of the desired sampling rate. As will be described in greater 
detail below, determining the derivatives of the function, in an accurate but approximate manner, 
5 facilitates approximating the function to determine an accurate function value for a specified 
sampling point at the desired sampling rate. More specifically, the estimates of the derivative 
values may be used to determine coefficients of a cubic polynomial function that is used to 
describe the unknown signal. 

To reduce the overall number of multiplications and additions required per unit time, the 
10 LPFs are all implemented in poly-phase form. As indicated in Figure 4, two-fold interpolation 
filters implemented in poly-phase form consists of two sub-components, referred to as "poly- 
phase components", operating at the low sampling rate followed by a MUX that selects between 
the poly-phase component outputs. The MUX operates at the high sampling rate. Figure 5 shows 
an example implementation of LPF 3 in poly-phase form. The outputs R0 and Rl correspond to 
15 poly-phase components 1 and 2, respectively of Figure 4. 

Figure 6 shows the top level block diagram of the CCSI block; the intermediate signal, f, 
with sampling rate 352.8 kHz enters the block and it generates output samples at a rate of 48 
kHz. 

Figure 7 is a functional block diagram that illustrates a clamped cubic spline interpolator 
20 formed according one embodiment of the present invention. The clamped cubic spline 
interpolator (CCSI) of Figure 7 includes a "Signal Buffer" block, an accumulator block, "ACC", 
a "Derivative Estimation" block, a "CCS Coefficient Calculation" block, and a "CCS 
Evaluation" block. The CCSI of Figure 7 is coupled to receive an upsampled signal that is to be 
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decimated by a specified non-integer amount. In the specific embodiment shown, the CCSI is 
coupled to receive a 352.8 kHz signal and to produce a decimated signal having a sampling rate 
of 48.0 kHz. More specifically, the signal buffer receives the input signal of 352.8 kHz and 
captures a plurality of samples which it produces in parallel. More specifically, the signal buffer 

5 captures 10 sequential samples of the 352.8 kHz input signal in the described embodiment of the 
invention. The signal buffer produces the 10 buffered samples to the derivative estimation block. 
The signal buffer also produces two select samples to the CCS coefficient calculation block. The 
derivative estimation block produces two derivative estimates based upon the 10 buffered 
samples received from the signal buffer. The two derivative estimates produced by the 

10 derivative estimation block are produced to the CCS coefficient calculation block. The CCS 
coefficient calculation block then produces, in the described embodiment, four coefficient values 
for a cubic polynomial. The four coefficient values are produced to the CCS evaluation block, 
which is further coupled to receive a value (x) corresponding to a specified decimated sampling 
point. The CCS evaluation block then produces a function value at the desired sampling rate of 

15 48 kHz. Operation and theory of each of the blocks within Figure 7 are described in greater 
detail below. 

A clamped cubic spline is a cubic polynomial of the form 



f(x) = a 0 + a x x + a 2 x 2 + a 3 x 3 , ( 1 ) 

20 

where the coefficients ao,..,<23 depend upon known values of the function and its 
derivatives in the boundary points jco and X\, and where the interpolation point x (the point at 
which an estimate of the underlying function is desired) is in-between the boundary points. For 
convenience,/^) is defined such thatx e [ 0 ; 1]. Specifically, 

25 
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«o =/(*o) 
«■ = /'(*o) 

a 2 =3{f(x,)-f(x 0 ))-{f'(x x ) + 2f'(x 0 )) 
a 3 =2(/(x 0 )-/(x l )) + /'(x 1 ) + /^o) 



(2) 



In the present application the signal is generally unknown and hence the derivatives of 
the signal are not known analytically a-priori. Thus, for high performance applications, precise 
5 estimates of the signal derivatives in the boundary points must be generated. The present 
invention uses Taylor approximation and Richardson extrapolation to arrive at highly precise 
derivative estimates. Details of the approach are given in a subsequent section. 

For discussion of Figure 7, assume that the block labeled "Derivative Estimation" in 
Figure 7 generates precise estimates of the signal derivatives in the boundary points in-between 
10 the appropriate samples of the intermediate signal. Specifically, these derivatives are evaluated 
as / \x 0 ) and / '(*/), and are labeled "fpO" and "fpl", respectively. The "CCS Coefficient 
Calculation" block then calculates the CCSI coefficients, as defined by (1) and (2). The 
evaluation of the polynomial corresponding to (1) is performed by the "CCS Evaluation" block. 
The block labeled "ACC" provides sampling triggers at an average rate of 48 kHz (the desired 
15 sampling rate in the described embodiment). 

A Richardson extrapolation is performed by the derivative estimation block for producing 
precise derivative estimates. Richardson extrapolation originates from the field of numerical 
analysis, where it is used to generate increasingly precise estimates of a function A(h) which 
depends upon a parameter h. For example, h could be the step size of a method for numerical 
20 integration or differentiation. In the following, it will be assumed that 

A(h)->A for A->0. 
Furthermore, with the assumption that A(h) is a "smooth" function of h, i.e., sufficiently 
many derivatives of A(h) exist for A=0, Taylor expansion yields, 

25 




(3) 



If k denotes the order of the method or function A{h), by definition, 
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fl , =a 2 =••• = <*,_, =0, (4) 

and thus 

A(h)=A + a k h lc +a m h m +o(h m + l ), (5) 
5 where the possibility 

a k + \ ="k + 2 = - = ««-i =0 (6) 
has been accommodated for. 

Now let hi and A 2 be two given A- values — specifically, for the purposes of this 
application, 

10 /z, = 2h 2 . (7) 

From (3), 

A(h i )=A + a k h l k + a m h; + 0(hr) 



A(h 2 )=A + a k h 2 k +a m h 2 m + 0(hr l ) 



It follows that the a*-term can be eliminated by multiplying the second equation with 
1 5 (h\/h 2 ) k = a k and subtracting: 

A _ a 'Ay-A(h,) + a _ a^_ h . +0{ltr) (9) 



or, since a=2, the order Richardson extrapolation is given by 



20 



A(H^A(k 2 ) A^-f'\ (10) 



2*-l 



Note that the fc* order Richardson extrapolation results in an order approximation of 
A. Thus, Richardson extrapolation may be repeated, this time to //1 th order, for even higher order 
25 approximations of A . 
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Consider next how Richardson extrapolation may be employed for highly precise 
derivative estimates in the evaluation of the clamped cubic spline. Forward Taylor 
approximation of a function^*) around the point x; yields 

f(x, +h)= f( Xi ) + h f'( Xi ) + \h 2 /'(*,)+ 1/. 3 /"(c), (1 1) 

where 

ce[x,. \x i +/?]. 

Backward Taylor approximation of a function^x) around the point x,- yields 

f(x i -h)=f(x i )-hf(x i ) + ih 2 r(x i )-ih 3 r(c), (12) 

where 

c g[x, -h \ x,]. 

Subtracting (12) from (11) eliminates even order derivatives, i.e., 

f( Xi + h)-f( Xi -h)= 2 h f'( Xl ,) + i/*7 m (<r), (13) 

where 

£e[x,-A ;x,. +/*]. 

It follows that 




For fixed x,-, introducing the compact notation 

f( Xi +jh) = f., j = ... -2,-1,0,1,2,... (15) 

and performing second order Richardson extrapolation of (14) according to (10) yields 
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/ >_ /-2-8/-,+8/ 1 - / 2 

Analogously, 



//= / - i " 8/ ; 2 + / /2 " /3 + 4 4 )- d7) 



Performing Richardson extrapolation of (16) and (17), this time to fourth order, yields 



f . -/-, + 40/_ 2 -256/,, +256/, -40/,+/, ( t) 

and 



- A + 40/,, - 256/ 0 + 256/ 2 - 40/ 3 + / s / 6 x 
1 360/* V 7 ' V ^ 

respectively. The 0(h 6 ) term is a remainder or error term. Because h has a value that is less 
than one, the term vanishes since h is being raised to the sixth power. In a practical sense, and as 
will be demonstrated shortly, the error is sufficiently small for hi-fidelity applications. 



It follows from equations (18) and (19) that calculating precise derivative estimates can 
be viewed as a filtering task. Specifically, (18) and (19) show that the task may be performed by 
an anti-symmetric FIR filter with impulse response 

h i z ) = m (l-40z' 2 +256z" 3 -256z" 5 +40z -6 -z -8 ). (20) 
To explain the above process in simple terms, the above approach attempts to provide 

reasonably accurate estimates of the derivatives by looking at surrounding information. The 

exponents of z reflect delay (-6 means six delay elements ago in the signal buffer of Figure 8, 

for example). The constant factors 1, 40 and 256 are mere multiplication terms. In the described 
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embodiments, as shown in Figure 9, only two sets of 3 (ml, m2 and m3) multiplication blocks 
are needed for implementing the algorithm of equation (20). 

Referring to Figure 8, a block diagram of the "Signal Buffer" appropriate for calculating 
derivative estimates using fourth order Richardson extrapolation is shown according to one 
embodiment of the present invention. Notice from equations (18), (19), and (20) a total of 10 
signal values are needed to estimate the derivatives f'o and /'/. Thus, as shown in Figure 8, the 
"Signal Buffer" block consists of 9 registers to hold the samples Jj, 

Table 2 (below) shows algorithm details of the "ACC" block of Figure 7. This block is 
driven by the 352.8 kHz sampling clock of the intermediate signal. The block maintains a 
modulo- 147 accumulator, which is incremented by 20 during each cycle. Generally, an overflow 
occurs whenever addition of the value of 20 is added to an amount within the accumulator which 
is initialized to zero but maintains a remainder from a previous overflow condition. For 
example, in a first pass, the accumulator is set to "0". After the 7 th addition of the value "20" the 
total value is still only 140. After the 8 th addition of the value "20", however, the accumulator 
value is 160 which is an overflow condition with a remainder of "13" which is equal to 160-147. 
For the next pass, the accumulator starts with the overflow remainder, 13, and adds the value of 
"20" thereto as the process continues. 

It can be shown that overflow of this accumulator indicates that a new output sample 
should be generated. Overflow occurs at an average rate of 48 kHz, corresponding to the desired 
output sampling rate. The signal "trig" is asserted when overflow occurs, indicating that the 
various processing blocks of Figure 7 must complete a processing cycle. The interpolation point, 
i.e., the point in the interval between xo and xj where the signal value is to be estimated, is 
found via a 20-entry look-up table (LUT). 
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ACC Operation 

set acc = 0; 

lutPointer = 0; 
loop 

acc = (acc + 20) modulo 147; 
trig = carry; 
if trig 

lutPointer = (lutPointer + 1) modulo 20; 
x = LUT(lutPointer); 
end 
end 

Table 2 - Operation of ACC ~ ~~ 

Figure 9 is a functional block diagram of the derivative estimation block according one 
5 embodiment of the invention. The symmetry of the Richardson extrapolation indicated by (20) 
can be exploited to reduce the number of multiplies to 3 rather than 6. These multiplies are 
shown as mi, m2 5 and m3 and are found by rewriting (20) as 

*M= M^-^hM^ -z- 2 )+m(z- 3 -z- 5 )- (2D 

10 

Generally, the derivative estimation block of Figure 9 calculates derivatives in the end 
points. The operation of Figure 9 is an implementation of formulas 18 and 19. As may further 
be seen, only three multipliers are necessary, namely, 1, 40, and 256. The outputs are estimates 
of the signal at endpoints xO and xl . 
15 The circuit of Figure 10 produces estimates of derivatives in the end points (the two 

points that straddle a sampling point at the desired (new) sampling rate. 

Figure 10 is a functional block diagram of the CCS coefficient calculation block that 
shows the details of the calculation of the CCS coefficients, according one embodiment of the 
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invention. The arithmetic indicated in this figure directly follows from equation (2). The factors 
of 2 and 3 are implemented as a left-shift and a left-shift and add, respectively, and are thus not 
considered multiplications in the described embodiment. The CCS coefficient calculation block 
receives the two buffered values fo and fj as well as the derivative estimate values fpo and f p i. 
5 Figure 1 1 shows the details of the CCS evaluation block of Figure 7. As may be seen, the 

CCS evaluation block receives the four polynomial values ao-a 3 as well a (x) from the 
Accumulator as described herein. 

The method of polynomial evaluation adopted here is referred to as Horner's method. 
This method is arrived at by noting that the cubic polynomial of equation (1) may be re-written 
10 as 

f(x) = (((a 3 x) +a 2 )x + a, )x + a 0 . (22) 
The advantages of Horner's method are a reduced number of multiplications necessary to 
evaluate the polynomial and typically more well balanced intermediate results. Thus, the circuit 
of Figure 1 1 implements equation (22) and produces a function output (magnitude of the random 
15 signal) at a time value of (x) which is the sampling point according to the desired sampling rate. 
Stated differently, the CCS evaluation block evaluates the function and produces a result at the 
desired sampling rate sampling times. 

Figures 12A and 12B are flow charts of a method of the present invention. More 
specifically, Figure 12A is a flow chart that illustrates a method for near unity sampling 
20 according to one embodiment of the present invention. Initially, the inventive embodiment 
includes receiving and up-sampling digital data representing a continuous waveform signal, the 
digital data having a first sampling rate (step 500) and filtering the upsampled data (step 502). 
The inventive method further and optionally includes receiving and up-sampling the upsampled 
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data to produce twice upsampled data and filtering the twice upsampled data (step 504). In the 
described embodiments, the invention further optionally includes receiving and up-sampling by a 
factor of two the twice upsampled data to produce thrice upsampled data and filtering the thrice 
upsampled data (step 506). As an alternative to upsampling by a total factor of eight times over 
5 three upsampling steps, one may upsample by a factor of 8 (or to the preferred upsampled 
amount) in one step. One advantage of utilizing multiple steps, however, relates to filter design 
and complexity to block images created from each increase in sampling rate. Larger upsample 
step sizes require sharper filters as is known by one of average skill in the art. Finally, the 
inventive embodiment of Figure 12A includes producing estimated derivative values, producing 

10 coefficient values and producing data with a decimated data sampling rate from a clampled cubic 
spline (CCS) interpolator to produce data at the second sampling rate wherein a ratio defined by 
the first and second data sampling rates is near unity (step 508). 

Figure 12B is a flow chart that illustrates a method for up-sampling digital data 
representing a continuous waveform signal, the digital data having an initial sampling rate, 

15 wherein the up-sampler produces upsampled and filtered data an upsampled rate (step 510). 
Thereafter, the inventive method includes receiving data at the upsampled rate, and determining 
a value of (x), the sampling point at the desired sampling rate, and further producing a trigger at 
the sampling points at the desired sampling rate (step 512). Additionally, the embodiment of the 
inventive method further includes receiving the upsampled and filtered data at the upsampled 

20 rate, the signal buffer for producing buffered samples that is N samples wide wherein N is a 
value determined by a Richardson formulation for generating derivative estimates (step 514). 
The previous step facilitates producing at least two derivative estimates based upon the buffered 
samples (step 516) and therefore producing polynomial coefficient values based upon the at least 
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two derivative estimates and upon at least two select buffered samples produced by the signal 
buffer (step 518). 

Figure 13 shows the magnitude response in dB of the sampling rate converter in the 
frequency range 0-30 kHz. It demonstrates the flatness of the passband as well as the steepness 
5 of the transition band. As stated before, the passband corresponds to the frequency range 0-20 
kHz, whereas the transition band is the frequency range 20 kHz-24 kHz. Figure 14 shows a 
close-up of Figure 13. 

Figure 15 shows the typical ripple of the passband magnitude response. Notice that this 
plot is a function of the passband magnitude response of all three lowpass filters, i.e., LPFj, 
10 LPF 2 , and LPF 3 . 

Figure 16 shows a typical simulation result for one embodiment of the present invention. 
Specifically, Figure 16 shows the output power spectral density (PSD) in units of decibel relative 
to the carrier (dBc) for a 15 kHz sinusoidal input signal The simulated signal-to-noise-and- 
distortion ratio (SNDR) is 105 dB and the spurious-free dynamic range (SFDR) is better-than 
15 108 dB. 

Figure 17 shows a block diagram of the scheme employed to evaluate the SNDR and 
SFDR. The digital input sinusoid of known frequency enters the sampling rate converter, here 
labeled x[n]. The converter output, labeled y[n], is analyzed by first subtracting an estimate of 
the input signal x[n] and then performing a fast fourier transform (FFT) to arrive at the power 
20 spectral density of the error signal. This process continues in a feedback loop until the error 
signal power is minimized, resulting in an optimal estimate of the signal x[n]. 

Figure 18 shows simulated SNDR and SFDR vs. frequency of the sinusoidal input signal 
x[n]. Better-than 96 dB SNDR and 100 dB SFDR is observed over the audio frequency band. 
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Figure 19 shows the simulated magnitude response of the sampling rate converter found 
as the ratio of output signal amplitude and input signal amplitude. It is in agreement with the 
frequency response plot of Figure 15. 

While the invention is susceptible to various modifications and alternative forms, specific 
5 embodiments thereof have been shown by way of example in the drawings and detailed 
description. It should be understood, however, that the drawings and detailed description thereto 
are not intended to limit the invention to the particular form disclosed, but, on the contrary, the 
invention is to cover all modifications, equivalents and alternatives falling within the spirit and 
scope of the present invention as defined by the claims. As may be seen, the described 
10 embodiments may be modified in many different ways without departing from the scope or 
teachings of the invention. 
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